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One of the most intriguing phenomena in na- 
ture is the fractional quantum Hall effect (FQHE) 
observed in the half-filled second Landau level [1] 
which, arising in even-denominator filling factors, 
v = 5/2 and 7/2, is completely different from other 
FQHEs in its origin, all of which, except for those 
two filling factors, occur in odd-denominator frac- 
tions. Usually formulated in terms of a trial 
wave function called the Moore-Read Pfaffian 
wave function {H, current leading theories 0-0] 
attribute the origin of the 5/2 FQHE to the for- 
mation of Cooper pairs, not of electron, but of 
the true quasi-particle of the system known as 
composite fermion [8]. The nature of supercon- 
ductivity resulting from such Cooper pairing is 
particularly puzzling in the sense that it appar- 
ently coexists with strong magnetic fields, which 
poses an interesting dilemma since the Meissner 
effect is the most important defining property of 
superconductivity. This apparent dilemma is re- 
solved by the fact that composite fermions do 
not respond to external magnetic field at even- 
denominator filling factors. To provide direct ev- 
idence that it is composite fermions that actually 
form the superconducting condensate, here, we 
develop a numerically exact method of creating 
a Cooper pair of composite fermions and explic- 
itly compute the superconducting order param- 
eter as a function of real space coordinates. As 
results, in addition to direct evidence for super- 
conductivity, we obtain quantitative predictions 
for superconducting coherence length. Obtain- 
ing such theoretical predictions can serve as an 
important step toward fault-tolerant topological 
quantum computation (91. 1 1 Of] . 

We begin by emphasizing that no assumption is made 
throughout this work regarding the nature of the 5/2 
ground state. Our point of view is that, if it exists, 
the superconducting order parameter should arise nat- 
urally when computed via exact diagonalization. Com- 
puting the superconducting order parameter, however, 
explicitly depends on our choice of the quasi-particle in 
the ground state. In the usual Fermi liquid, the quasi- 
particle is identified with something very close to the free 



electron so that the superconducting order parameter is 
computed in terms of the "free" electron operator. The 
situation is less clear for strongly correlated electron sys- 
tems. It is, thus, useful to provide independent, unbi- 
ased evidence for superconductivity that is not suscepti- 
ble to our knowledge of the quasi-particle nature. Such 
evidence is obtained in the ground state energy itself. 
Since a prerequisite for superconductivity is the pairing 
between quasi-particles, the ground state energy must os- 
cillate depending on whether the particle number is even 
or odd. 

In order to obtain the FQHE ground state energy, it 
is convenient to use a compact geometry with periodic 
boundary condition. In this study, we use the so-called 
spherical geometry where electrons are confined on the 
surface of a sp here which has a magnetic monopole at its 
center ll|, Il2j . The total magnetic flux piercing through 
the spherical surface is given by 2Q4> with fa, called the 
flux quantum, being equal to 2irhc/e and with Q, called 
the monopole strength, being equal to either an integer or 
a half integer. Angular momentum eigenstates in this ge- 
ometry are known as the monopole harmonics which form 
a basis set to span the single-particle Hilbert space of each 
Landau level. Many-particle FQHE ground states are 
obtained by numerically diagonalizing the inter-particle 
interaction matrix within specific Landau levels; for ex- 
ample, the diagonalization is done within the half-filled 
second Landau level for the 5/2 FQHE. It is important 
to note that, due to the nature of the monopole harmon- 
ics, the ratio between the number of particles, N, and 
that of flux quanta, 2Q, is not exactly equal to the ther- 
modynamic filling factor, v, in finite-size systems; N/2Q 
approaches v as N —¥ 00. It is known from various stud- 
ies that the 5/2 FQHE state occurs at 2Q = 2N - 3 with 
an apparent "flux shift" of —3 

Figure [T] shows the results of the ground state energy 
per particle for the Coulomb interaction, He , of the half- 
filled second-Landau-level (SLL) in comparison with that 
of the half-filled lowest-Landau-level (LLL) as well as 
that of a model two-body interaction, H2, for which the 
Moore-Read Pfaffian wave function is known to be very 
close to the exact ground state [141 ] . As one can see, 
even-odd effects are clearly visible for both SLL He and 
H2 (top panels of Fig. [I]), while the same is not true for 
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FIG. 1: Even-odd effect The ground state energies per par- 
ticle are computed for the PfafEan two-body interaction, H2, 
[top left panel], the Coulomb interaction, He, projected in the 
second-Landau-level (SLL) [top right panel] and the lowest- 
Landau- level (LLL) [bottom two panels]. For the Coulomb 
interaction, the energy is given in units of e 2 /el oc (where 

loo — Ib\J2vQ/N with Ib = y^hc/eB []]|) and the usual 
charge background correction is taken into account. In the 
case of the LLL He, two different sectors of flux quantum 
number, 2Q, are studied for a given particle number, N: the 
Pfaffian flux of 2Q = 2N — 3 and the composite-fermion-sea 
flux of 2Q = 2N - 2. 



LLL Hq (bottom panels). In fact, the behavior of the 
LLL He ground state energy is quite different from the 
SLL counterpart, showing a peculiar pattern of repeating 
local minima at N = 4, 9, and 16 (bottom right panel 
of Fig. Q]) . The energy minima would repeatedly appear 
whenever N becomes n 2 with n being a positive integer if 
numerical diagonalization is possible for bigger systems. 

This behavior can be understood as follows. It is well 
established by now that the ground state of the half- 
filled lowest Landau level is the composite fermion (CF) 
sea which, while being a compressible state in the ther- 
modynamic limit, has a shell structure in terms of the 
CF Landau levels in the finite-size spherical geometry. 
The ground state energy exhibits local minima whenever 
energy shells below the chemical potential become com- 
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with n being a 



Note that the correct flux sec- 
tor for the CF sea is obtained at 2Q = 2N — 2 for reasons 
provided by the CF theory. One may wonder what hap- 
pens if the same flux sector studied for SLL Hq, i. e., 
2Q = 2N — 3, is investigated in the LLL. The LLL re- 
sults are plotted in the bottom left panel of Fig.[TJ where 
it is shown that energy minimum positions are a little bit 
different from those of the CF sea flux sector, while still 
not exhibiting the even-odd effect. Here, energy min- 
ima occur for a combination of two reasons. Now that 
the flux is shifted from that of the CF sea, composite 
fermions are subject to small, residual field, resulting in 
two ways of minimizing the total energy. The first is 
to minimize the kinetic energy of composite fermions by 



filling the similar shell structure as before, but with a 
different particle number sequence of N = n(n + 1) with 
n being a positive integer. The second is to reduce any 
residual interaction energy by forming Wigner crystals 
of composite fermions, which occur at N = 4, 6, 8, 12, 
and 20 with N corresponding to the number of vertices 
in regular polyhedrons of Euclid. The energy reduction 
due to the second effect is seen via a weak local mini- 
mum at N = 8 and a kink at N = 4. The conclusion is 
that the SLL He definitely exhibits the even-odd effect, 
which strongly supports the existence of superconductiv- 
ity in a complete analogy with what happens in small 
superconducting grains [17| . On the other hand, there is 
no such evidence for the lowest Landau level in any flux 
sectors. The clear signature of the even-odd effect in the 
half-filled SLL, but not in the LLL, is one of the most 
important results in this work. It is worthwhile to men- 
tion that the difference between the second and lowest 
Landau-level physics stems from very subtle quantitative 
changes in the Coulomb matrix elements parameterized 
by the Haldane pseudopotential. 

Emboldened by strong support from the even-odd ef- 
fect, we now set out to address the main issue in this 
paper: what are the quasi-particles that are being paired 
and what would be the appropriate order parameter for 
them? At this stage, it is important to distinguish be- 
tween two types of the quasi-particles in the system. The 
first is the aforementioned composite fermions that, as 
shown in the following, form Cooper pairs in a similar 
fashion as electrons in the BCS theory with exactly the 
same statistics. The second is vortices that, in a direct 
analogy to their counterpart in the p + i p super conduc- 
tivity, may possess non-Abelian statistics [18j, [l_9( . In this 
study, we are interested in the superconducting nature of 
composite fermions. 

The superconducting order parameter is convention- 
ally defined in the BCS theory as the expectation value 
of two fermion creation operators, say, and c_ k , in the 
fixed-phase coherent ground state. For fixed-number sys- 
tems, an alternative but physically identical formulation 
of the BCS superconducting order parameter is given in 
terms of the matrix element between the ground states 
in the N and N + 2 systems [loj]. In this work, we use 
such formulation and compute the superconducting or- 
der parameter in a numerically exact manner without 
any (BCS or otherwise) approximation. Instead of the 
usual momentum space representation, it is much more 
convenient here to compute the superconducting order 
parameter in real space: F(r) = (N + 2|c^(r)ct(0)|A r ) 
where |iV + 2) and \N) are the ground states in the iV + 2 
and N systems, respectively. The preceding definition, 
however, is not yet appropriate for the 5/2 FQHE since 
the true elementary quasi-particle is not an electron, but 
a composite fermion [21[ . A necessary condition for creat- 
ing a composite fermion is that the magnetic flux should 
be increased by 24>q, i. e., Q — ¥ Q + 1, at the same time as 
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an electron is created. For this reason, a correct defini- 
tion for the 5/2 FQHE superconducting order parameter 
may be given as follows: 

F*(r) = (yV + 2,Q + 2|4(r)4(0)|iV,Q) , (1) 

where \N + 2, Q + 2) and \N, Q) are the ground states 
for the particle-flux sector of (N + 2, Q + 2) and (N, Q), 
respectively. In the above, ct indicates the composite 
fermion creation operator accompanied with an appro- 
priate increase in flux. We emphasize that Eq. ([1]) is con- 
sidered an appropriate order parameter since it is propor- 
tional to the superconducting gap within the BCS the- 
ory [20| . It vanishes in the normal state not by symmetry 
reasons, but due to gap vanishing. Finally, it is interest- 
ing to note that the composite fermion pair operator used 
in Eq. (Q} is conceptually similar to the composite boson 
operator considered by Read [22j and by Rezayi and Hal- 
dane [23|, which exhibits Bose-Einstein condensation at 
odd-denominator fractions. 

The composite fermion creation, however, requires 
more than just an increase in total flux: two flux quanta 
should be locally bound to an electron. In fact, the rea- 
son for the apparent absence of the Meissner effect is that 
composite fermions do not respond to external magnetic 
fields at half filling since all magnetic fields are captured 
away in the very process of forming composite fermions. 
Recalling that a flux quantum is actually equivalent to 
a correlation hole with 2tt phase twist, what one has to 
do is to make sure that appropriate correlation holes are 
bound to each added electron. Unfortunately, however, it 
is close to impossible to derive the CF creation operator 
as a function of that of electrons in a closed analytic form 
since a composite fermion is an extended object contain- 
ing full many-body correlations. In this study, we take 
a different approach. Since we do not know a priori how 
correlation holes are attached to the constituent electron 
inside a composite fermion, we let the system evolve to 
reach the lowest energy configuration. In other words, we 
perform exact diagonalization in the increased particle- 
flux sector of (N + 2, Q + 2) under the constraint that 
an electron is pinned at the origin and the other at the 
position of r. Let us call the lowest energy state in this 
situation |?/;^ r+2,(3+2 (r, 0)). Relegating the discussion for 
details to the following paragraph, we show the density 
profile of +2 < Q+2 (r, 0)) in Fig.H 

Let us now discuss how to obtain +2 < Q+2 (r, 0)). 
To this end, it is necessary to know the LLL-projected 
operator for creating an electron at the position of r, 
which is given by: 




where r = i?(sin 9 cos (j>, sin 9 sin <j>, cos 9 — 1) with the ra- 
dius R = IbVQ an( i the magnetic length lg — y/Tic/eB. 
YQi m is a complex conjugate of the monopole harmonics 
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FIG. 2: Creation of a composite-fermion Cooper pair 

a, Schematic diagram for the finite-size spherical geometry 
used in our numerical computation, b, Resulting density pro- 
file of a composite-fermion Cooper pair as a function of sep- 
aration distance in the system of N — 14 (which includes 
two electrons participating in the created Cooper pair). Den- 
sity profiles are computed in both cases of H2 (top panels) 
and the SLL He (bottom panels) for comparison. The far 
left panels depict the situation where two composite fermions 
(bright yellow spots) are most separated, i. e., one located in 
the north pole and the other in the south pole. On the other 
hand, the far right panels show that two composite fermions 
are almost on top of each other. Note that, with separation 
distance being measured along the chord, figures show the 
correct distance while the shape is distorted. 

with the monopole strength, Q, the total angular momen- 
tum, I, and the z-componcnt angular momentum, m. 
is the usual operator for creating an electron in the m- 
eigenstate. The LLL-projected density operator is then 
defined by p pro j(r) — c Pro j( r ) c proj{r) with normalization 
that the maximum value of density is unity. Equipped 
with the LLL-projcctcd density operator, one can ob- 
tain |</>f +2 ' Q+2 (r,0)) by using the Lagrange multiplier 
method. That is, the constraint of pinning electrons at 
the origin and the r position are imposed by introducing 
the Lagrange multiplier term in the Hamiltonian: 

H = H int - M[p P roj{r) - 1] - A 2 [/5 prO j(0) - 1] , (3) 

where Hi nt is the inter-particle interaction Hamiltonian. 
The Lagrange multipliers, Ai and A2, are determined 
via imposing the saddle-point conditions for the ground 
state energy: de gr /dXi = de gr /d\2 = 0. In the end, 
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|'0^ r+2 ' < ^ +2 (r, 0)) is given as the ground state satisfying 
the above saddle-point conditions. It is worthwhile to 
mention that, when two electrons are separated across 
the north and south poles, the constraint can be im- 
posed precisely without recourse to the Lagrange mul- 
tiplier method. We have verified that the results from 
the Lagrange multiplier method are completely consis- 
tent with that of the precise method. 

The final ingredient for computing the superconduct- 
ing order parameter of the 5/2 FQHE is the realization 
that |f/;^ r+2 ' < ^ +2 (r, 0)) obtained from exact diagonaliza- 
tion is actually a normalized wave function: 



N+2,Q+2 



(r,0)> 



4(r)4(0)|iV, Q) 



yj (TV, Q|c,(0)c*(r)ct(r)ct (0)| JV, Q) 



where normalization is insured by the denominator which 
is nothing but the square root of the pair distribution 
function for composite fermion. Arguing that the CF pair 
distribution function should be similar to the electronic 
counterpart, g(r), we arrive at the following conclusion: 
cUr)ci(0)\N,Q) = y^)|Vf +2 ' Q+2 (r,0)). The CF su- 
perconducting order parameter is finally written as fol- 
lows: 



F*(r) = y/g(rj(N + 2,Q + 2\1> 



N+2.Q+2 



(r,0)>. (5) 



Figure [3] shows F* (r) for both H 2 and the SLL He ob- 
tained from exact diagonalization of TV = 8 — 18 systems, 
which are compared with that of the LLL He- As one 
can see, for Hi, the superconducting order parameter ex- 
hibits a rather smooth curve, being robust across various 
finite-size systems while, for the SLL He, there are some- 
what more fluctuations even though the overall behavior 
is similar. The superconducting order parameter is essen- 
tially negligible for the LLL He, which, along with the 
absence of the even-odd effect, is consistent with the re- 



cent conclusion of Storni, Morf, and Das Sarma |24j. It is 
emphasized that all our numerical results taken together 
convincingly establish that the real 5/2 FQHE state has 
underlying superconductivity. In addition to providing 
concrete theoretical evidence for superconductivity, an 
important lessen from F*(r) is that the superconduct- 
ing order parameter is rather long-ranged, being sizable 
at least up to ten magnetic lengths. This is important 
since the superconducting coherence length sets a natural 
length scale for performing coherent quantum operations. 

Now, let us make a connection between our findings 
and previous results from field theories. The compos- 
ite fermion (CF) creation operator defined in our paper 
creates a bound state between an electron and a corre- 
lation hole under the condition that the magnetic field 
is increased by two flux quanta whenever an additional 
electron is added. In the lowest Landau level, a corre- 
lation hole is always concomitant with the formation of 
a vortex. Since the magnetic field is increased by two 
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FIG. 3: Superconducting order parameter as a func- 
tion of separation distance between two constituent 
composite fermions in a CF Cooper pair The top and 

bottom panels show the results for H2 and the SLL He, re- 
spectively. Here, TV indicates the total number of electrons 
including the two participating in the created Cooper pair. 
Note that, for the TV = 16 and 18 systems, the Hilbert space 
becomes too large to diagonalize except when two added com- 
posite fermions are separated across the north and south 
poles. Dashed lines is a guide to eye. In order to accentu- 
ate the distinction between the superconducting FQHE state 
at v — 5/2 and the usual non-superconducting states, we 
compute the superconducting order parameter for the ground 
state of the half-filled lowest Landau level, the results of which 
are denoted by down triangles in the figure. Note that only 
two data points, TV = 6 and 18, are available for the half- 
filled lowest Landau level since the CF sea state is unsta- 
ble unless the number of electrons comprising the CF sea 
(which does not include two electrons participating in the 
added Cooper pair) is equal to that of the filled shell struc- 
ture, i. e., TV — 2 = n 2 with n being a positive integer. Also, 
only even particle numbers are relevant for possible pairing. 



flux quanta, the created bound state is actually equiv- 
alent to that between an electron and a double vortex. 
Meanwhile, in the composite fermion Chern-Simons (CS) 
gauge field theory, it is generally accepted that binding 
a double vortex is adiabatically connected to the attach- 
ment of an infinitesimally thin tube of two flux quanta, 
i. e., Chern-Simons flux, into an electron. In this sense, 
our CF creation operator essentially plays the same role 
as creating CS-flux attached electrons. According to 
the theory of Halperin, Lee, and Read [25|, and also 
Kalmeyer and Zhang [26j , the ground state at v = 1/2 is 
the Fermi sea state of the CS-flux attached electrons and 
therefore, in this state, the superconducting order param- 
eter must vanish when measured in terms of the CS-flux 
attached electrons. On the other hand, the Moore-Read 
Pfaffian state, which is the exact ground state of H3, by 
construction contains the superconducting order in terms 
of the CS-flux attached electrons. It is for the first time 
in this work that these expectations are explicitly shown 
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to be true in realistic Hamiltonians. 
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